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ABSTRACT 

The measurement of radio polarization may be regarded as a three- 
dimensional statistical problem because the large photon densities at radio wave- 
lengths allow the simultaneous detection of the three Stokes parameters which 
completely describe the radiation's polarization. The statistical nature of the 
problem arises from the fluctuating instrumental noise, and possibly from fluc- 
tuations in the radiation's polarization. A statistical model is used to derive 
the general, three-dimensional statistics that govern radio polarization measure- 
ments. The statistics are derived for specific cases of source-intrinsic polarization, 
with an emphasis on the orthogonal polarization modes in pulsar radio emission. 
The statistics are similar to those commonly found in other fields of the physi- 
cal, biological, and Earth sciences. Given the highly variable linear and circular 
polarization of pulsar radio emission, an understanding of the three-dimensional 
statistics may be an essential prequisite to a thorough interpretation of pulsar 
polarization data. 

Subject headings: Methods: statistical - polarization - pulsars: general - radio 
continuum: general 



1. INTRODUCTION 



The large photon densities at radio wavelengths, unlike the much smaller densities at 
much shorter wavelengths, allow the simultaneous detection of the three Stokes parameters 
which completely describe the polarization of the electromagnetic radiation (Radhakrishnan 
1999; Thompson, Moran, & Swenson 2001). The randomly fluctuating instrumental noise in- 
troduces a statistical element to the measurements, so a complete and proper interpretation 
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of radio polarization observations requires a three-dimensional statistical analysis. Interest- 
ingly, but quite understandably, radio astronomers do not fully exploit this unique aspect of 
radio science, and effectively reduce a polarization measurement to a problem of one or two 
dimensions by investigating the circular or linear polarization, respectively, on an individual 
or separate basis. 

The polarization signature of the radiation mechanism under investigation usually dic- 
tates which Stokes parameters are actually measured. For example, the Zeeman effect pro- 
duces a distinctive, circularly-polarized signal, and the measurement of a single component 
of the polarization vector, the Stokes parameter V, is usually sufficient to observe the po- 
larization of the effect. Synchrotron radiation produces significant linear polarization, but 
no circular polarization (e.g. Rybicki & Lightman 1979), so most polarization observa- 
tions of synchrotron-emitting sources only utilize the Stokes parameters Q and U. It follows 
that the measurement of circular polarization is a one-dimensional (ID) statistical problem 
because only one Stokes parameter is involved. Similarly, the measurement of linear po- 
larization is a two-dimensional (2D) statistical problem because two Stokes parameters are 
used. The statistics applicable to the measurement of circular and linear polarization are 
well-understood, and are documented in McKinnon (2002). 

Three-dimensional (3D) statistics of radio polarimetry are particularly relevant to pulsar 
radio emission, where both the linear and circular polarization can be significant and highly 
variable. The radio emission mechanism of pulsars is not completely understood, and single 
pulse polarization observations are made in an attempt to understand it. The results of these 
observations are usually displayed as longitude-dependent histograms of fractional linear 
polarization, fractional circular polarization, and position angle (e.g. Backer & Rankin 
1980; Stinebring et al. 1984) or as joint probability densities of linear polarization and 
position angle (PA) or of circular polarization and PA (e.g. Cordes, Rankin, & Backer 
1978). The histograms and probability densities beautifully illustrate the statistical nature 
of the emission's polarization fluctuations and the fascinating orthogonal polarization modes 
(OPM). However, these display methods, by their very nature, force one to interpret the 
emission's linear and circular polarization on a separate basis (e.g McKinnon 2002), and 
reveal little, if anything, about the association of the two polarizations. The display methods, 
in effect, betray our reliance upon the familiar, but limited, ID and 2D polarization statistics. 
Applying the more general, and more descriptive, 3D statistics to pulsars may be essential 
to understanding their emission mechanism. One reason why the 3D statistics may not have 
been applied is it appears they have yet to be determined. 

Directional statistics are a subset of the more general 2D and 3D vector statistics, and 
govern the possible orientations of a vector. They are particularly applicable to experiments 
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where a vector's orientation, but not its amplitude, can be measured, such as measurements 
of the vanishing angles of birds taking flight and measurements of magnetic remanence in 
rock formations (e.g. Mardia 1972). The statistics' widespread applicability in biology and 
the Earth sciences, the pioneering statistical work of R. A. Fisher, and the advent of the 
modern computer have all contributed to the rapid development of directional statistics over 
the last 30 years (Fisher, Lewis, & Embleton 1987). Interestingly, directional statistics were 
originally developed to solve a problem in astronomy. In 1734, Bernoulli used the statistics 
to prove that the orbital planes of the planets in the Solar System could not be aligned by 
chance (Mardia 1972, and references therein). The statistics are not widely used in radio 
astronomy now, but the possibility of their application to radio polarimetry is very intriguing. 

The purpose of this paper is to derive and summarize the 3D statistics that are applicable 
to radio polarimetry. In §2, the statistics of the measured polarization vector are derived 
for a source of fixed polarization using the statistical model described in McKinnon (2002). 
The specific cases of no polarization, linear polarization only, and circular polarization only 
are evaluated. The polarization vector statistics for OPM are derived in §3. The directional 
statistics of a polarization vector are derived in §4 for the cases of fixed polarization and 
OPM. The applicability and limitations of the statistics are discussed in §5. The 3D statistics 
will be applied to polarization observations in a future paper. 



2. FIXED POLARIZATION 

The polarization of any source of electromagnetic radiation is completely described by 
the Stokes parameters Q, U, and V. The parameter V describes the circular polarization 
of the radiation, and the parameters Q and U describe the real and imaginary parts, re- 
spectively, of its linear polarization. The Stokes parameters form a right hand, Cartesian 
coordinate system, and the radiation's polarization can be represented by a vector in this 
three dimensional space, which shall be designated as "Poincare space" in what follows. For 
the general case of a polarization vector with an amplitude p and an orientation specified by 
an azimuth <p and a colatitude 9 Q , the circular polarization is fiy = pcos9 Q , and the real and 
imaginary parts of the linear polarization are hq = psin# o cos0 o and fijj = psm9 sm(f) , 
respectively. 

The Stokes parameters that are actually recorded during an observation include the 
contributions from the Gaussian instrumental noise of the radio telescope. If Xn is a ran- 
dom variable that represents the instrumental noise, the measured Stokes parameters are 
(McKinnon 2002) 
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Q = ^ Q + X NiQ , (1) 
U = fiu + X NiU , (2) 

V = /iy + X N>V . (3) 

Thus denned, each of the measured Stokes parameters is a Gaussian random variable with 
mean \x and standard deviation uat, where cat is the magnitude of the instrumental noise. 
Equations 1 through 3 are a set of linear equations that form a simple, statistical model 
for the polarization measured with a radio telescope. In its simplest form, the model com- 
pletely ignores any systematic effects such as instrumental imperfections, calibration errors, 
parallactic angle rotation, interstellar scintillation, and Faraday rotation. 

As long as the antenna temperature produced by the radio source is small in comparison 
to the system temperature of the radio telescope, the instrumental noise in each of the 
measured Stokes parameters is independent from that in the others. This being the case, 
the measured Stokes parameters are independent of one another, and their joint probability 
density is the product of their individual probability densities (i.e. Jquv — fQfufv)- The 
Stokes parameters can be evaluated in the context of a joint probability density because they 
can be detected simultaneously at radio wavelengths and thus the polarization state of the 
radiation can be completely determined at any instant. 

The possible amplitudes and orientations of the measured polarization vector can be de- 
termined by converting the joint probability density from Cartesian to spherical coordinates. 
If the random variables x, y, and z represent the measured Stokes parameters Q, U, and 
V, respectively, the measured amplitude of the polarization vector is r = (x 2 + y 2 + z 2 ) 1 ^ 2 , 
and its measured orientation 2 is given by a colatitude of 9 = arccos(z/r) and an azimuth 
of = arctan(y/:r). The position angle (PA) of the linear polarization is ip = 0/2. After 
converting from Cartesian to spherical coordinates, one can show that the joint probability 
density governing the measured values of the polarization vector's amplitude and orientation 
is 
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2 This nomenclature for orientation angles is consistent with that of traditional, spherical coordinate 
systems. It is different from the nomenclature used in McKinnon & Stinebring (1998) and McKinnon (2002) 
where azimuth and PA were denoted by 9 and 0, respectively. 
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2.1. No Polarization 

If the electromagnetic radiation is not polarized, the amplitude of its intrinsic polar- 
ization vector is p = 0. From equation 4, the joint probability density for the measured 
polarization vector is 

The distribution of the vector's amplitude is found by integrating the joint probability density 
over all orientation angles. 

f(r) = f /(r,MW= i\fl e M-&) <® 

JO JO a N V n V Z<J N/ 

The amplitude distribution is identical in functional form to the Maxwell-Boltzmann distri- 
bution for the speed of molecules in an ideal, classical gas. In retrospect, this is not surprsing 
considering that similar assumptions, namely the statistical independence of the Stokes pa- 
rameters and of the vector components of a molecule's velocity, were made in the derivation. 
The mean, (r), and standard deviation, a, of the amplitude distribution are 

(r) = 4<r N /Vto, (7) 

a = a N ^3 - 8/tt. (8) 
The distribution of the vector's orientation is 

f(0,<j>) = J™ f(r,9,<l>)dr = ^mne, (9) 

which is isotropic. 



2.2. Linear Polarization 

When the radiation is linearly polarized, the colatitude of the intrinsic polarization 
vector is o = tt/2. If the vector azimuth is <p = 0, the polarization occurs in the Stokes 
parameter Q, and the intrinsic Stokes parameters are fiQ = p and nu = nv — 0. The 
joint probability density for the measured values of this polarization vector's amplitude and 
orientation is 



sin^ 



T% (27T)3/2 



cxp 



(r 2 + p 2 ) 



exp 



rp sin 9 cos <fi 



'N 



(10) 



- 6- 



The distribution of the vector's amplitude is 
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The amplitude distribution is shown in the top panel of Figure 1 for different signal-to-noise 
ratios (s = p/<Jn)- As can be seen in the figure, the amplitude distribution is Maxwell- 
Boltzmann for small s, and resembles a Gaussian when s is large. The mean amplitude 
is 



(r) = jf rf(r)dr = a N ^-exp^- 
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In the limiting cases of high and low signal-to-noise ratio, s, the mean amplitudes are 

(r) ~p + a 2 N /p, p»a N , 



(12) 
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p < a N . 
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2.3. Circular Polarization 

The colatitude of the radiation's polarization vector is o = when it is circularly 
polarized. The Stokes parameters of the radiation's polarization are nv — P and /j,q = 
= 0. The joint probability density for the measured values of this polarization vector's 
amplitude and orientation is 
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(15) 



By integrating over all orientation angles, one can show that the distribution of the vector's 
amplitude is identical to the result in equation 11. Likewise, the mean amplitude of the 
measured vector is given by equation 12. It follows that one cannot distinguish between 
linear and circular polarization on the basis of the amplitude distribution alone. 

The joint probability density given by equation 15 is independent of and symmetric in 
azimuth. Therefore, the orientation of the polarization vector varies primarily in colatitude, 
9, and is distributed according to 

"2tt poo 

(16) 
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Colatitude, 6 (degrees) 

Fig. 1. — Distribution of a polarization vector's measured amplitude (top panel) and co- 
latitude (bottom panel) for radiation with fixed circular polarization. The distributions are 
shown for different signal-to-noise ratios, s = p/cr N . The colatitude distribution evolves from 
an isotropic distribution at s = to a Rayleigh distribution for large s. The value of the 
instrumental noise used in the amplitude distributions is = 1- 
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The colatitude distribution has only one free parameter, s, and is shown in the bottom panel 
of Figure 1. When s is large (s >> 1), the colatitude distribution approaches the Rayleigh 
distribution. 



Orthogonal polarization modes are commonly observed in the radio emission from pul- 
sars (e.g. Manchester, Taylor, & Huegunin 1975; Cordes, Rankin, & Backer 1978; Backer & 
Rankin 1980; Stinebring et al. 1984). In a rather confusing use of the term, the fact that the 
modes are orthogonal actually means their polarization vectors are antiparallel in Poincare 
space. The modes also appear to occur at the same time (McKinnon & Stinebring 1998, 
hereafter MSI; McKinnon & Stinebring 2000, hereafter MS2). So although the individual 
modes may be highly polarized, their simultaneous interaction can cause the observed polar- 
ization to be small. To account for the statistical fluctuations observed in the polarization of 
pulsar radio emission, the amplitudes of the mode polarization vectors can be represented by 
the random variables Xi and X 2 (MSI; MS2; McKinnon 2002). The means of these random 
variables shall be designated as ji\ and fi 2 , and their standard deviations shall be denoted 
as o\ and a 2 . Since the observations indicate that the modes occur simultaneously and their 
polarization vectors are antiparallel, the intrinsic polarization fluctuations occur along a di- 
agonal in Poincare space. For a diagonal specified by an azimuth <p and a colatitude 6 a , the 
observed Stokes parameters can be described by (MS2; McKinnon 2002) 




(18) 



3. 



ORTHOGONAL MODES OF POLARIZATION 



Q = sin# o cos0 o (Xi - X 2 ) + X NiQ , 



(19) 




(20) 
(21) 



V = cos^(X 1 -X 2 )+X NiV . 



The equations presented so far are valid regardless of the type of random variable that 
actually describes the amplitudes of the mode polarization vectors. The probability density of 
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Fig. 2. — Distributions of a polarization vector's measured amplitude (top panel) and colat- 
itude (bottom panel) when circularly-polarized OPM occurs in the radiation and the mean 
intrinsic polarization is p = 0. In this case, the radiation's polarization fluctuates, and the 
different curves represent the distributions expected for different mode fluctuation ratios, p. 
The value of the instrumental noise used in the amplitude distributions is = 1. 
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these random variables must be known before the joint probability density can be derived. In 
the examples which follow, the mode polarization amplitudes are assumed to be independent, 
Gaussian random variables. 



3.1. Linearly-Polarized Modes 

Consider the case when the modes are linearly polarized (9 = n/2) and their polariza- 
tion fluctuations occur in Stokes Q (<f> = 0). In this case, the fluctuations in the measured 
values of U and V are due to the instrumental noise. Therefore, they are Gaussian ran- 
dom variables with zero means (/iu — /ly — 0) and standard deviations of ay = ay = a N . 
Since the mode polarization amplitudes are assumed to be Gaussian random variables, the 
measured Stokes Q is also a Gaussian random variable with mean /j,q — p — \i\ — /j,2 and a 
standard deviation of oq = (of + of + ojv) 1 / 2 = <7jv(l + p 2 ) 1 ^ 2 - The quantity p is a measure 
of the intrinsic polarization fluctuations relative to the instrumental noise. With these as- 
sumptions, the modes will occur at PAs of ip = 0° and ip = 90°, and will occur with equal 
frequency if the mode polarizations are the same on average (i.e. /ii — = 0). The 

measured Stokes parameters are again independent of one another, and their joint probabil- 
ity density is the product of their individual probability densities (see, also, Appendix A3 of 
MSI). When the modes have identical mean values (/ii = /z 2 ; Pq = 0), the joint probability 
density is 



t 2 sin i 

/(r '"' 0) = ^(2^(iT7)^ exp 



r 2 (l + p 2 ) — r 2 p 2 cos 2 sin 2 9 



2(1 + P 2 K 

and the distribution of the measured polarization vector's orientation is 



(22) 



sin 9 1 + p 2 

4tt [1 +p 2 (l -cos 2 0sin 2 #)] 3 / 2 ' 



f{9A)= I f(r,9,<j>)dr = — H ■ — Zf rT - T — . (23) 



3.2. Circularly-Polarized Modes 

Although the observations indicate that OPM is predominantly linearly polarized, con- 
sider an idealized case where OPM is circularly polarized. The assumption of circularly 
polarized modes, combined with the observational facts that the modes occur simultane- 
ously and their polarization vectors are antiparallel, means that the intrinsic polarization 
fluctuations occur only in the Stokes parameter V (9 = 0). Again assuming that the mode 
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polarizations are Gaussian random variables with identical mean values, the joint probability 
density of the measured polarization vector's amplitude and orientation is 



fir A' 
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(24) 



The assumption of circularly polarized modes introduces azimuthal symmetry to the 
3D statistics, and allows for a straightforward, analytical derivation of the amplitude distri- 
bution. 



/(r)= SvSS exp (-^)i' 

The amplitude distribution is shown in the top panel of Figure 2 for different values of p. 
The figure shows that the modes tend to stretch what was originally the Maxwell-Boltzmann 
distribution of the instrumental noise to a highly asymmetric distribution. The extent of the 
tail on this distribution increases with the magnitude of the mode fluctuations, p. 

The azimuthal symmetry of equation 24 also means the orientation of the polarization 
vector varies primarily in colatitude. 



p 2 r 2 u 2 
2(1 + P 2 K, 



du 



(25) 



2 (1 + p A sin 6) 6 ' z 

The colatitude distribution is shown in the bottom panel of Figure 2 for different values of p. 
The distribution is isotropic for p = 0, and the polarization modes become more obvious as 
p increases. The two peaks in the colatitude distribution do not occur precisely at 9 = 0, 7r 
because the instrumental noise in Stokes Q and U weights the orientation of the polarization 
vector towards 9 = ir/2. 

Equations 23 and 26 were derived for polarization fluctuations along specific diagonals 
in Poincare space. For fluctuations along any other diagonal specified by a colatitude 9 and 
an azimuth O , the possible orientations of the measured polarization vector are described 
by 

f(a a\ sin g 1 + p 2 

/(m) = 7 mi + p mm)] 3 / 2 ' (27) 

where 



g(9, = 1 — [cos^ 2 cos 2 9 + sin 2 O sin 2 ^(cos 2 cos 2 O + sin 2 0sin 2 O )]. 



(28) 
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4. DIRECTIONAL STATISTICS 

The directional statistics of a polarization vector predict the possible orientations of the 
vector at a given vector amplitude, r a . They are determined from the 2D or 3D statistics 
by computing the conditional density of the vector's orientation. In the specific 2D case 
when only the linear polarization is considered, the conditional density gives the possible PA 
values on a circle of radius r . In the more general 3D case, the conditional density gives the 
possible vector orientations on a sphere of radius r Q . 



4.1. Distribution on a Circle 

Figure 3 shows a linear polarization vector in the Q-U plane of Poincare space, and 
illustrates the calculation of the PA conditional density. The measured vector is the vector 
sum of the radiation's polarization and the randomly polarized instrumental noise. The 
conditional density is the distribution of the measured vector's PAs along the circle of radius, 
r Q , and is given by 



U\r a (lp) 



f{ro,i/>) 



(29) 



For the geometry displayed in the figure and noting that ip is defined on the interval — 7r/2 < 
tp < 7r/2, the joint probability density and the amplitude distribution are, respectively 
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(30) 



(31) 



where Io(x) is the modified Bessel function of order zero. From Equations 29 through 31, 
one can show that the PA conditional density is a von Mises distribution (Mardia 1972; see 
also MSI) with a dispersion parameter of b = r Q fiQ/a 



N ' 



U\M) = 



exp(6cos2^) 



-tt/2 < i) < ir/2. 



(32) 



The PA conditional density is shown in the top panel of Figure 4 for different values of b. It is 
a uniform distribution at b = 0, and is a Gaussian distribution with a FWHM of (2 In 2/6) 1 / 2 
when b is large . 
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Fig. 3. — Graphical representation of the calculation of the PA conditional density. The 
radiation's polarization vector has an amplitude /j,q and points in the direction of increas- 
ing Q in the Q-U plane of Poincare space. The measured polarization vector includes the 
instrumental noise, which is a vector with an amplitude characterized by cn that can point 
in any direction (dashed circle). The PA conditional density is the distribution of measured 
PAs along the circle of radius r Q . 
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Fig. 4. — Conditional density for the position angle (PA) of a linear polarization vector. 
When the radiation's polarization is fixed, the conditional density follows the von Mises 
distribution (top panel). For linearly-polarized OPM, the conditional density is a bimodal 
von Mises distribution (bottom panel). The distributions are shown for different dispersion 
parameters, b. 
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For linearly-polarized OPM, the joint probability density is (MSI) 
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where p = (af + a^Y^/a^. When the modes have the same polarization on average (pq = 0), 
the amplitude distribution is 
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One can then show that the PA conditional density for linearly-polarized OPM is a bimodal 
von Mises distribution with a dispersion parameter of b = r 2 p 2 /Aa 2 N (l + p 2 ). 

exp[6 cos(4'0)] 



7T/ (6) 



-tt/2 < i) < ir/2 



(35) 



The bimodal von Mises distribution is shown in the bottom panel of Figure 4 for different 
values of p. For display purposes, equation 35 was modified by replacing ip with ip — n/4 so 
that the modes occur at if) — —it/A, 7r/4 in the figure. 

Detection thresholds are commonly applied to the measured amplitude of the linear 
polarization vector before constructing PA histograms (e.g. Stinebring et al. 1984; MSI). 
The resulting histogram is not an exact replica of the underlying PA distribution, but is 
more similar to a PA conditional density. 



4.2. Distribution on a Sphere 

When all three Stokes parameters are measured, the conditional density of a polarization 
vector's orientation is the distribution of vector orientations on the surface of a sphere, instead 
of the perimeter of a circle as in § 4.1. As with the PA conditional density, the conditional 
density of the vector's orientation is calculated from the joint probability density and the 
amplitude distribution. 

For a radiation source with fixed polarization properties, the joint probability density is 
given by equation 4 and the amplitude distribution is given by equation 11. Since the distri- 
bution of the polarization vector's amplitude is independent of its orientation as discussed 
in § 2, one can use the two equations to show that the conditional density always follows the 
general Fisher distribution, 

f(6,(f)\r ) = ^( r °' ^lA. — KsinO exp{/4cosfl cosfl + sin6> o sin6>cos(0 - (f> )]}, (36) 
j(r ) 47rsmh(/c) 
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where k = pr /a 2 N is a concentration parameter (Fisher, Lewis, & Embleton 1992). As a 
specific example, 9 Q = when the radiation is circularly polarized, and from Equation 36 
the conditional density is 

f(9, (f)\r Q ) = KSm6 exp(re cos 9). (37) 
An smh(/t) 

Equation 37 is independent of azimuth, as one would expect from the azimuthal symmetry 
of the problem. Therefore, integration of f(9, 4>\r ) over yields a factor of 2tt, and the 
resulting colatitude conditional density is fo\ ro {9) = 2nf(9 : (j)\r ). The colatitude conditional 
density is shown in the top panel of Figure 5 for different values of k. Equation 37 can 
also be derived directly from the joint probability density in Equation 15 and the amplitude 
distribution in Equation 11. 

For the specific case of circularly-polarized OPM, the conditional density can be found 
from Equations 24 and 25. The conditional density is a Watson bipolar distribution (Fisher, 
Lewis, & Embleton 1992) with a concentration parameter of k — r 2 p 2 /2(l + p 2 )u%, 

f(9,<P\r ) = -^-eM^os 2 9). (38) 

Here, w(k) is a normalization factor given by 

w(k) = / exp(KU 2 )du. (39) 
Jo 

This particular distribution is also independent of azimuth, so the colatitude conditional 
density is again fe\r a {9) = 2irf(9,4>\r ). The colatitude conditional density is shown in the 
bottom panel of Figure 5 for different concentration parameters. Both the Fisher and Watson 
bipolar distributions become isotropic distributions in the limit k — 0. 



5. DISCUSSION 

The essential point of this paper is microwave polarization measurements can be in- 
terpreted in the context of a joint probability density of the Stokes parameters Q, U, and 
V because the large photon densities at radio wavelengths allow their simultaneous detec- 
tion. These same photon statistics make coherent phase detection possible and thus the 
measurement of visibility phase for radio aperture synthesis observations (Radhakrishnan 
1999). This important realization allows the derivation of the general 3D statistics for a 
measured polarization vector. The 3D statistics, in turn, can be used to extract informa- 
tion about the radiation's polarization properties. The statistics derived here are unique to 




50 100 150 



Colatitude, 9 (degrees) 



Fig. 5. — Conditional density of a polarization vector's colatitude, 6, when the radiation 
is circularly polarized. When the radiation's polarization is fixed, the conditional density 
follows the Fisher distribution (top panel). In the case of circularly-polarized OPM, the 
conditional density follows a Watson bipolar distribution (bottom panel). The distributions 
are shown for different concentration parameters, k. 
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radio polarization measurements, particularly when the radiation's polarization fluctuates 
on a timescale comparable to the sampling interval. The statistics may not apply to polar- 
ization observations at much shorter wavelengths because the Stokes parameters cannot be 
measured simultaneously and the polarization state of the radiation may change between 
consecutive measurements. 

The 3D statistics provide the appropriate framework for interpreting observations where 
both the linear and circular polarization are significant and/or highly variable, as is often the 
case in pulsar radio emission. Historically, but with some notable exceptions (e.g. Manch- 
ester, Taylor, & Huguenin 1975), the linear and circular polarization of individual pulses 
have been interpreted or evaluated separately, making it difficult to visualize or determine 
how the polarization vector's orientation varies from pulse-to-pulse. By applying the 3D 
statistics derived in this paper, one may be able to make more definite statements about the 
true polarization state of the radiation. 

The statistical model used to derive the polarization statistics accurately describes most 
radio polarization observations. The model is a set of linear equations that accounts for the 
radiation's polarization properties and includes the instrumental noise. The model is valid 
as long as the telescope system temperature, Tg, is greater than the antenna temperature 
produced by the radio source, Ta- When Ta becomes comparable with Tg, the instrumental 
noise is not independent between Stokes parameters and the model is no longer linear. 

In deriving the 3D statistics, special emphasis was placed on the specific case of circularly- 
polarized radiation. This apparent bias does not signify a shortcoming of the statistical 
model's ability to describe other types of polarization, but merely highlights the fact that 
straightforward, analytical solutions can be derived given the azimuthal symmetry provided 
by the geometry (e.g. Eqns. 17 and 38). Also, the functional form of the distribution 
of the polarization vector's orientation is easier to illustrate (e.g. Figs. 1, 2, and 5) for 
circularly-polarized radiation because the distribution contains only one independent vari- 
able (colatitude). 

Both the 3D and directional statistics for microwave polarimetry are remarkably similar 
to those that are relevant in other fields of the physical, biological, and Earth sciences, and 
it is astonishing that such diverse measurements can be described by the same statistics. For 
example, measurements of remanent magnetization in rock formations are consistent with 
the Fisher distribution, and measurements of the intersections between the cleavage and 
bedding planes of certain types of minerals follow the Watson bipolar distribution (Fisher, 
Lewis, & Embleton 1987, and references therein). The von Mises distribution can describe 
the orientations of sand-grains and the vanishing angles of birds taking flight (Mardia 1972, 
and references therein). These similarities hold great promise for future polarization work 
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because the mature, rigorous, statistical methods developed in other fields may be applied 
to radio polarization observations. 
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